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Abstract. We present recent investigations on dipole and quadrupole excitations in spherical skin 
nuclei, particular exploring their connection to the thickness of the neutron skin. Our theoretical 
method relies on density functional theory, which provides us with a proper link between nuclear 
many-body theory of the nuclear ground state and its phenomenological description. For the calcu- 
lation of the nuclear excited states we apply QPM theory. A new quadrupole mode related to Pygmy 
Quadrupole Resonance in tin isotopes is suggested. 
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INTRODUCTION AND THEORETICAL MODEL 

The development of the new experimental facilities on radioactive beams moved forward 
the nuclear structure physics to the regions of nuclei far from stability yj]. One of the 
most interesting results was the discovery of a new dipole mode at low-excitation energy 
in nuclei with neutron skin [Q]]. This bunching of 1 states resembling spectral structures 
known to indicate resonance phenomena was named Pygmy Dipole Resonance (PDR). 
The presence of a neutron skin should affect also low-energy quadrupole states. Due to 
oscillations of neutrons from the surface region one may expect a formation of a soft 
quadrupole mode, namely Pygmy Quadrupole Resonance (PQR), which resembles the 
properties of the PDR. 

Here, we present our investigations on the dipole and quadrupole excitations in nuclei 
from N=50, 82 and Z=50 regions. The approach is based on a Hartree-Fock-Bogoljubov 
(HFB) description of the ground state B2|] by using a phenomenological energy-density 
functional (EDF) [Hi 3]. The excited states are calculated with the Quasiparticle-Phonon 
model (QPM) pD. The model Hamiltonian is given by: 

H = H MF + Hj* + Hf M + (1) 

which has the structure of the traditional QPM model [HI] but in detail differs in the 
physical content as discussed in [4]. The HFB term Hmf = H sp + H pa j r contains two 
parts: H sp describes the motion of protons and neutrons in a static, spherically-symmetric 
mean-field. H pa i r accounts for the monopole pairing between isospin identical particles 
with phenomenologically adjusted coupling constants, fitted to data Q6Q. The mean- 
field Hamiltonian compromises microscopic HFB effects [0, Hi and phenomenological 
aspects like experimental separation energies and, when available, also charge and mass 




FIGURE 1. (left) Proton (dash line) and neutron (solid line) ground states densities in Sn isotopes; 
(right) One-phonon QPM calculations for the total PDR strengths in the 100 ~ 132 Sn isotopes (upper panel) 
are displayed for comparison together with the nuclear skin thickness 8r (lower panel). Experimental data 
on the total PDR strengths in 116 Sn and 124 Sn 10] and 112 Sn JUl are also shown. In the lower panel, the 
skin thickness derived from charge exchange reactions by Krasznahorkay et al. ll9l [Toll 



radii. For numerical convenience the mean-field potential is parametrized in terms of 
(superpositions of) Wood-Saxon (WS) potentials with adjustable parameters [0]. 

The remaining three terms present the residual interaction H res = Hj^ + H^ M + H^f . 
As typical for the QPM, we use separable multipole-multipole H^ 1 and spin-multipole 
Hg M interactions both of isoscalar and isovector type in the particle-hole and multipole 
pairing Hfif in the particle-particle channels, respectively. The isoscalar and isovector 
coupling constants are obtained by a fit to energies and transition rates of low-lying 
vibrational states and high-lying GDR's fflfllll. 

The nuclear excited states are described by quasiparticle-random-phase- approxima- 
tion (QRPA) phonons. They are defined as a linear combination of two-quasiparticle 
creation and annihilation operators as follows: 

qU = \l (vjfAinUf) - tifanU/j) , (2) 

where j = (nl jmx) is a single-particle proton or neutron state; At and are time- 
forward and time-backward operators, coupling proton and neutron two-quasiparticle 
creation or annihilation operators to a total angular momentum X with projection /i by 

means of the Clebsch-Gordan coefficients C^., m , = (jmfm'\Xfi) . 

The phonons obey the QRPA equation of motion, which solves the eigenvalue prob- 
lem, i.e. giving the excitation energies E{ and the time-forward and time-backward am- 
plitudes [5] wj*i 2 and <p^ 2 , respectively. 



E*= 0- 8 MeV E"= 8- 8.5 MeV . „ E*= 8- 20 MeV 
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FIGURE 2. (left) QRPA proton (dash line) and neutron (solid line) dipole transition densities summed 
over the 1 ~ excited states in a given energy region ( indicated in the figure ) in neutron-rich Sn isotopes 
and (right) in 88 Sr and 90 Zr, respectively. 



As discussed in [0, LL2l LL3] we obtain important additional information from the 
analysis of the spatial structure of the nuclear response on an external electromagnetic 
field. That information is accessible by considering the one-body transition densities 
8p(r), which are the non-diagonal elements of the nuclear one-body density matrix. 
The transition densities are obtained by the matrix elements between the ground state 
|*F ? -) = \JiMi) and the excited states l^f) = \JfMf). The analytical procedure for the 
calculation of the transition densities is presented in details in Q4£| . 



RESULTS 

We have pointed out in our previous papers [0, 0] that the proper determination of the 
ground state of the nucleus is essential for the extrapolations of the QPM calculations 
into unknown mass regions. Here, we put special emphasis on a reliable description 
of the mean-field part, reproducing as close as possible the g.s. properties of nuclei 
like binding energies, root mean square radii, and separation energies to the required 
accuracy along an isotopic chain. This is achieved by solving the ground state problem 
in a semi-microscopic approach IHIU]. The calculated ground state neutron and proton 
densities are displayed in Fig.l (left) for several tin isotopes. Of special importance for 
our investigation are the surface regions, where the formation of a skin takes place. The 
skin thickness is defined with the deference between neutron and proton rms radii with 



We identify in QPM |/jM/) = |0) with phonon vacuum and obtain the excited states by means of the 
QRPA state operator, eq. ©, \J/M f ) = Q\ .|0). 




FIGURE 3. (left) QRPA calculations of quadrupole states in Sn; (right) Proton (dash line) and 
neutron (solid line) quadrupole transition densities summed over the QRPA 2 + excited states in a given 
energy region ( indicated in the figure ) in l20 Sn. 



the equation 8r = y < r 2 > — < r 2 > and its evolution with the mass number A for 

the tin izotopic chain is presented in the lower panel of Fig. 1 (right) (see also in [0]). 
From both figures (Fig. 1 (left) and Fig. 1 (right) lower panel) follow that for A> 106 the 
neutron distribution begins to extend beyond the proton density and the effect continues 
to increase with the neutron excess, up to l32 Sn. Thus, these nuclei have a neutron skin. 
The situation reverses in 100_102 Sn, where a tiny proton skin appears. In investigated 
N=50, 82 isotones (N>Z) the presence of neutron skin was observed as well | |14Lll5|1 . 

Of s pec ial interest are the low-energy dipole excitations in tin [3, 4] and N=50, 82 
nuclei [|l4l . \l5 ]. which we studied in respect to their connection with skin vibrations. 
From QRPA calculations in 112 ~ 132 Sn and N=50, 82 nuclei we obtain a sequence of low- 
lying dipole states located below the neutron threshold of almost pure neutron structure. 
These states exhaust only a tiny fraction, less than 1% of the total Thomas -Reiche-Kuhn 
(TRK) energy weighted dipole sum rule and are related to PDR [0, SLED, In general 
the PDR states are of mixed symmetry and their total transition strength is closely related 
to the size of the neutron skin [4] . The dependence of the calculated total PDR strength 
on the mass number in 100 ~ 132 Sn nuclei is presented in comparison to 5r in Fig.l (right). 
The total PDR strength increases when 8r increases and correspondingly the neutron or 
proton skin thicknesses increase. 

The analysis of the dipole transition densities acts as further indicator of PDR and 
provides a precise separation of the PDR states from the GDR. The calculated QRPA 
proton and neutron dipole transition densities for several Sn isotopes and N=50 iso- 
tones (N>Z) are presented in Fig. 2. At excitation energies below the neutron emission 
threshold ( presented in the first column of the left and right plots of Fig. 2 ) we observe 
in-phase oscillation of protons and neutrons in the nuclear interior, while at the surface 
only neutrons contribute. This picture is associated with PDR mode. The states above the 
neutron threshold (second and third column of the left and right plots) indicate isovector 



oscillations corresponding to the excitation of the GDR. An interesting finding is the 
most exotic 100 Sn nucleus where the lowest-lying 1 state has a proton structure. It has 
been related to a proton PDR and investigated in details by analysis of dipole transition 
densities |4|]. 

We perform QRPA calculations of 2 + states in 120 Sn up to excitation energies E*=35 
MeV. Similar to the PDR we observe a sequence of states, located just above the collec- 
tive 2|" and below the ISGQR (see Fig. 3 (left)) with predominantly neutron structure. 
The specific admixture of isoscalar and isovector components in their phonon content 
is a signature of mixed symmetry states (the amplitudes of proton and neutron state 
vectors are with opposite signs). They are connected with the ground state with small 
B(E2) transitions. By analysis of proton and neutron quadrupole transition densities (see 
Fig 3 (right) bed ) these states can be related to neutron PQR, where neutrons from 
the periphery region oscillate against isospin symmetric core. In addition we can dis- 
tinguish precisely the PQR states from the collective isoscalar 2^ (Fig. 3 (right)a) and 
ISGQR (Fig. 3 (right)e). In Fig. 3 (right)f QRPA neutron and proton transition densities 
corresponding to the IVGQR are presented. From the calculation the positions of the 
maximums for the ISGQR and IVGQR in 120 Sn are determined, which are E*= 12.2 
MeV and E*= 28.1 MeV, respectively. These values are close to the given ones in [fl6ll . 

In conclusion QPM calculations in Z=50 and N=50, 82 nuclei predict low-energy 
dipole strength in the energy region below the particle emission threshold directly 
related to the size of the neutron(proton) skin |4(]. The unique character of this mode 
is associated with PDR, further confirmed by the shape and structure of the transition 
densities clearly distinguishable from the GDR. The discussed results are in a good 
agreement with the experimental data [0, 03, H]. In addition quadrupole states were 
investigated in 120 Sn. By analyzing their structure and proton and neutron transition 
densities clear separation between collective isoscalar, isovector and mixed symmetry 
states is achieved. A new quadrupole mode related to PQR in 120 Sn, resembling the 
properties of PDR is suggested. Supported by DFG project Le 439/1-4 and BMBF 
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